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Abstract 

We present a model of iV-component hard-core bosons on a lattice. The limit N — > oo can be solved 
exactly. A saddle point approximation leads to a 1/iV expansion and allows the calculation of physical 
quantities like the density of the condensate, the correlation between the components and the density- 
density correlation function and the correlation between different components. We find a superfluid phase 
with a tendency towards a Mott insulator at high densities and finite N. 

1 Introduction 

A system of interacting lattice bosons undergoes a phase transition between a superfluid and a Mott 
insulator ^ |2 E] . Motivated by the recent discovery of a Mott-insulating state in an ultracold gas of 
Rubidium atoms in an optical lattice |JJ, we study the effect of strong (i.e. hard-core) interaction in a 
gas of lattice bosons. The strength of the interaction is tuned effectively by introducing N components 
in the Bose gas [5-12]. These components can be considered as N degenerate states per site of the optical 
lattice, assuming that there is a hard-core interaction only in the same state. Thus, the bosons can avoid 
hard-core interaction by choosing an unoccupied state, leading to an effective interaction between the 
bosons that becomes weaker with an increasing number of states N. Another tunable parameter is the 
tunneling rate t between neighboring sites of the optical lattice. For sufficiently small tunneling rate the 
Bose-Einstein condensate disappear, since for t = (atomic limit) there can be no condensate but only 
a (Mott) insulator. In order to study the properties of such a multicomponent Bose gas we evaluate 
physical quantities like the particle density, the condensate density, the structural order parameter with 
respect to the N components and the density-density correlation function. Our method is a saddle-point 
integration for large values of N and a 1 /N expansion ^3] • 

2 The Model 

The Hamiltonian for a Bose gas in a <i-dimensional optical lattice with Af sites and N states per site is 

N 

r,r' a,a' = l 

where a\ a (a r ^ a ) is the creation (annihilation) operator of a hard-core boson at site r on the optical 
lattice and in state a, i.e. (a£ Q ) 2 = 0. The states a = 1,2, TV refer to N different components of the 
Bose gas. t ra y a i is the tunneling rate for a hard-core boson from the state a' at site r' to the state a 
at site r. Here we will assume isotropic tunneling 



, _ f J/N for r, r' nearest neighbors 

L otherwise 
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For the functional integral we introduce 

(w r>r >) = e t/2 

with t = {t r ,a,r'a') (cf. |15p. 

This system can be characterized by a diagonal order parameter for the components a: 
which is the number operator of component a. An off-diagonal order parameter is 

a l,a a r,a> («' ^ a) (1) 

which measures how likely it is to switch a particle from component a' to component a. Besides the 
expectation value of these order parameters we also need their correlations to describe ordering of these 
degrees of freedom. 

2.1 Functional Integral and Saddle-Point Approximation 

The starting point is the grand-canonical ensemble of the iV-state hard-core Bose gas. It is defined by 
the partition function Zhc 

Z HC =Tre- fj{H -^ nr \ 

where —fj, ^ r n r controls the particle number in the grand-canonical ensemble 13 .In the following, we 
will use the fugacity C = gAI instead of the chemical potential \i. It was shown JTJ] that Zhc can 
be expressed by a functional integral with respect to two complex fields tp x and Xx ■ x — (r, t) are 
coordinates of a space-time lattice, where the time is t — 1, 2, ...,/?. The field is periodic with respect to 
time: ifp+i, r = The partition function then reads 

f i r d<p x dtp*dxxdx* x „_s / 9 n 
Zhc = znJ[[ ^~y 2 e (2) 



with the action 



S:=|£(l + C<->*^ + f El**! 2 



and the prefactor 



-EE'^ + ^+^fe+yi (3) 

x a. 

z N = {Aetv^)N 2pN '. 

The number of states N only appears in front of the TV-independent action in the partition function. 
Therefore the functional integration can be performed for large values of TV by a saddle-point (SP) 
approximation 15 j. The corresponding SP equations can be solved with uniform solutions: 

fx =<P, Xx= X- (4) 
With Y, x > Kx> = \ tw0 uniform SPs are found, a trivial SP 

V = = X = X* = (5) 

labeled by v — and a non-trivial SP 

M = 2^1^ ~- 

iX = ~f, iX* = -\v* (6) 

with broken U(l) symmetry, labeled by v = 1. The complex fields ip x and Xx are expected to fluctuate 
about the SP solution 

<p x = <p + 5ip x , ip* x = <p* + 5cp* x , Xx = X + $<f>x , X* x = X* + 5x1 ■ ( 7 ) 
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In order to study the effect of the fluctuations, we introduce a four component real field 

§ , , = 5<t>2,v _ 5<Px 
$fo,x ~ $X'x 

V 6^ x J \ 6-C J 
Substituting this result into the partition function leads to a Gaussian approximation 

Zhc,v ~ Z v Zqp^ u , 
where the first factor is the partition function for N ~ oo 

Zq = zn 



Z\ = zn exp 



-JVflA/^l-i-logCj 



(8) 

(9) 

(10) 
(11) 



and the second factor describes the partition function of quasiparticles, which arise from 1 /TV-corrections 



JQP,v 



ndS<f> x 



>) 



(12) 



where (...,...) denotes a scalar product. 

The prefactor f3J\f in Eq. comes from the summation over all lattice sites x. The Fourier 

components of the quasiparticle Green's functions are 4 x 4-matrices of the form 



-iA„ 



l + A v 



(13) 



with the 2 x 2-matrices 



vZ 1 . 



BI + 1 + 2B k cos(wr) 



B k cos(ujt) + 1 B k sin(wr) 
— B k sin(wr) B k cos(wr) + 1 



B k = 1 — r J + r-j Y^ d j=i cos kj an d 



Ao = C 



1 
1 



Ax = 



2/C-l 
1 



Now we can perform the Gaussian integral over the fluctuating field (|12(l and get 

Z QP ^=N-^\{det(gll 2 ) kuJ . 



(14) 



2.2 Spectrum of quasiparticles 

In the superfluid phase, the matrix {Gi)k ui can be recognized as Green's function of quasiparticles, the 
elementary excitations out of the condensate. The quasiparticle poles of the eigenvalues of det(^ 1 )fc w 
lead to the excitation spectrum. Identifying icu with the excitation energies e, one gets for small wave 
vectors k: 

" = ^«- 1 >" + (5* a ) > 

In the region (« 1, where /iq w 1 + ( is valid for the chemical potential, one can identify m = d/ J and 
gets the Bogoliubov spectrum for the weakly interacting Bose gas 




The linear dispersion relation for small wave vectors k is due to a Goldstone mode. 
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3 Physical quantities 



3.1 Correlation functions 

Correlation functions can be defined by introducing the generating fields 77" 1 and r]^ 2 and a generating 
functional 

where 

S g en({vf}, {vf} ■■= j + MtvWl' + J E l^l 2 

x,x' x 

-EE + + + < x )(^ + % + )] • 

a; a 

77" 1 generates the field of the bosons and ry" 2 the field of empty sites. According to Eqs. (0 and J2J the 
partition function is given as 

z HC = z gen ({vf},{vf})\^ nS ^ ■ 

An rt-particle correlation function can be defined as 

1 d 2n Z, 



C {X\Ot \ , . . . , x n (y. n 5 X n/?/! , . . . , X n$l) 



gen 



Zhc d^l 1 . . . drjg-drfe 2 . . . drff 



r;° 1 =?7j 2 =0 



All physical quantities that we present in the following are defined in terms of these expectation 
values. The only exception is the total density of particles, that can be defined as 

n := 1 — C(xa]xa). (15) 

All physical results are evaluated for zero temperature. 

3.2 Total density and density of the condensate 

The total particle density per site and component l|15fl gives in the large- N limit 

l_ C -i if C> i • (16) 

The 1 /iV-correction of n in the empty phase £ < 1 vanishes but in the superfluid phase the total density 
is enhanced. For large values of £ the density approaches the value 1 (the "full" lattice), indicating the 
tendency to a Mott insulating state. 

Our definition of the condensate density is based on off-diagonal long-range order, using the single 
particle correlation function 



n := lim - V C((r, 0)a; (r' , t)a) (17) 

r' — *r*r^ 1 * ■ 



r—r — >oo _ 

t=l 



^ /3 

For this quantity, the large- N limit is 

J if C<1 njn 

n ° - { (- 1 - r 2 if c > 1 ( J 

and indicates a depletion for large values of £. The 1/A-corrections even show a destruction of the 
condensate density for sufficiently large £, which can be regarded as another hint at the formation of a 
Mott insulator, where the phase coherence and therefore the condensate is expected to be destroyed (see 
figures 1 and 2). 

Because of the Goldstone mode in the excitation spectrum of quasiparticles, the 1/A-corrections 
diverge for the dimensions d = 1 and d = 2, in agreement with the expectation that in these cases there 
should not be a condensate. 
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3.3 Structural order parameter 

The appropriate quantity describing the correlation of the internal structure of the components of a given 
site is given by the expectation value of (JIJ as 

S a ,/3 := C(xa; x/3) . (19) 

For a = (3 this is equal to 1 — n, n being the total density of particles. For a ^ [3 the N — > oo-limit is 
equal to that of the condensate density, but the 1/iV-correction is slightly higher. As for the condensate 
density, the corrections indicate a destruction of the condensate for large values of £• Therefore this 
correlation function seems to be suppressed in a Mott insulating phase, too (see figure 3). 

3.4 Density-density correlation function 

The relation l|15|l suggests the definition of a density-density correlation function as 

V al3 {x, x) := 1 - C(xa; xa) - C(x x'0) + C(xa, x'0; x'/3, xa) . (20) 
In the large- TV limit this correlation function reduces to the square of the particle density: 

lim V af3 (x;x') = n 2 

The correlation length of spatially decaying 1/iV-corrections diverges at the critical fugacity £ = 1: 

6,«(c-i)-* • 

4 Conclusions 

We have found a depletion of the condensate caused by strong interaction in a dense ./V-component Bosc 
gas. This effect leads eventually to the destruction of the condensate and to the formation of a Mott 
insulator for sufficiently large fugacity £ at N < oo. However, the condensate survives in the limit 
N — > oo, since the interaction is always weak in this case. On the other hand, in the very dilute regime 
(i.e. for C ~ 1) the condensate density increases with decreasing N. This reflects the well-known fact that 
an increasing interaction supports BE condensation in a dilute Bose gas The structural order 

parameter S a ,/3 measures the correlation between the N different states of the Bose gas. This quantity 
behaves qualitatively like the condensate density but differently on a quantitative level: the dilute regime, 
where the interaction supports this parameter, is broader and its destruction happens at higher values of 

C- 

Our saddle-point approximation is similar to the Bogoliubov theory at low densities |16| but deviates 
substantially from the latter at higher densities. In particular, the appearence of a gap in the excitation 
spectrum of the Mott-insulating phase is beyond the Bogoliubov theory. 
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Appendix 

Evaluation of correlation functions 

Correlation functions can be written as expecation values with respect to the fields tp, x- The results for 
the functions mentioned in this work are: 



C{xa,x'(3-,x'(3,xa) = ([l + Q{ip x +i Xx ){w* x +ixl)] 1 [1 + C(<Px> + *X*')(¥>*' + ***')] ') , «) ¥= {?' , P) 



The functions can be expanded up to second order in the fields ip and Xi and the SP approximation is 
applied. The zeroth order provides the large- N limit and the second order leads to a 1/A-correction that 
appears as linear combination of two different expectation values: 



C{xa- xa) = ([1 + C(v* + iXx)(¥* x + iX* x )] *) 




,(x,a) ± (x',(3) 



$(r - r', t 



((S<p x + i8xx){5<Px> + iSXx')) = {{&<fil + ^X x ){5^ + iSx* x >)) 




e i[k(r-r')-u>(t-t')] 



i_ r d d k r 

N J J2^J 



dw B(k)(B(k)-(cos(uj)) 

2^ 2B(k) cos(w) - C + B(k) 2 (( - 2) 



e i[fe(r-r')-«(t-t')] 



with the function B(k) := 1 — J + ^ X^j=i cos 

Spacially constant corrections depend on $(0, 0) and $(0,0) only. The integral with respect to uj can 
be performed exactly and the fc-integration was done numerically. 
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Figure 1: Total particle density and condensate density for J = 0.1 
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Figure 2: Condensate fraction n /n for J = 0.1 
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Figure 3: Structural order parameter for J = 0.1 
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